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Abstract 

In general, quantum field theories (QFT) require regularizations and infinite renormalizations 
due to ultraviolet divergences in their loop calculations. Furthermore, perturbation series in the- 
ories like QED are not convergent series, but are asymptotic series. We apply neutrix calculus, 
developed in connection with asymptotic series and divergent integrals, to QFT, obtaining fi- 
nite renormalizations. While none of the physically measurable results in renormalizable QFT is 
changed, quantum gravity is rendered more manageable in the neutrix framework. 
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The procedure of regularization and renormalization is a big step forward in making sense 
of the infinities that one encounters in calculating perturbative series in quantum field the- 
ories. The result is a phenomenal success. For example, Quantum Electrodynamics (QED), 
the paradigm of relativistic quantum field theories, suitably regularized and renormalized, is 
arguably the most accurate theory ever devised by mankind. Yet in spite of the impressive 
phenomenological successes, the specter of infinite renormalizations has convinced many, 
including such eminent physicists as Dirac and Schwinger, that we should seek a better 
mathematical and/or physical foundation for quantum field theory, without simultaneously 
tearing down the towering edifice we have built on the existing one. In another development, 
Dyson has shown that the series as defined by the Feynman rules in QED is not a con- 
vergent series and has suggested that it is instead an asymptotic series in the fine structure 
constant a, i.e., in the number of internal integrals (for given outside lines). In this paper, 
we propose to apply neutrix calculus, in conjunction with Hadamard integrals, developed 
by J.G. van der Corput |2| in connection with asymptotic series and divergent integrals, to 
quantum field theories in general, and QED in particular, to obtain finite results for the 
coefficients in the perturbation series. (A more detailed discussion^ will appear elsewhere.) 
The replacement of regular integrals by Hadamard integrals in quantum field theory appears 
to make good mathematical sense, as van der Corput observed that Hadamard integrals are 
the proper tool to calculate the coefficients of an asymptotic series. (Actually Hadamard 
integrals work equally well for convergent series.) 

We begin by recalling the definition of asymptotic series 0]. The series f(x) = ao + ai(x — 
b) + a 2 (x — 6) 2 + ... for finite b is an asymptotic series if and only if there exists an n > 0, 
such that for n > n , 



with the remnant being at most ~ (x — b) n+1 . 

Next, following van der Corput |2], we define a neutrix as a class of negligible functions 
defined in a domain, which satisfy the following two conditions: (1) the neutrix is an additive 
group; (2) it does not contain any constant except 0. Let us illustrate the concept with the 
following example considered by Hadamard: for s real, 




(1) 
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For s > 0, the integral converges even as £ — > 0. For s < 0, Hadamard neglects £ s /s and 
log£ as £ — > 0. Here we have a neutrix which we will call N(0), consisting of functions 
KO = e (0 + c i£ s + c 2 l°g£) where e(£) — > as £ — > 0, and where C\, c 2 , and s are arbitrary- 
constants. This results in writing 

2 1 f s-^ 5 for s 4 

X s " 1 cfc = <^ (3) 
"(P) [log 2 fors = 

Note the analytic extension in the complex s plane of the answer for Re s > to the entire 
complex plane with the exclusion of s = 0. 

Before applying neutrices to QED, we need to consider the generalized Hadamard neutrix 
H a defined to contain the negligible functions 

v(0 = U(0 + e(0, (4) 

where e(£) — > as £ — > a. Each of the functions U(£) is defined by an asymptotic series 
based on a: 

U(0 ~ ^ =oXh (£ - af» log fc *(£ - a). (5) 

Here \hi and integers kh > are independent of £, Re ^ — >■ oo as h — > oo, and log fch (£) 
stands for (log(£)) fc, \ An example is provided by 

f {z - a)^ 1 log k (z -a)dz = {k + l)" 1 \og k+1 {b - a). (6) 

Similarly, we can define the Hadamard neutrix by Eq. (j^J) where now e(£) — > as £ — > oo 
and the function [/ (£) has a Hadamard development in powers of £ _1 in its asymptotic series: 

U(0 ~ ££° =oX/l log^ £, (7) 

where Re — >■ — oo as h — > oo. 

We can now demonstrate a very valuable property of the Hadamard neutrix. Recall 
that in the theory of distributions developed by Schwartz, generalized functions usually 
cannot be multiplied. Consider, for example, the one-dimensional Dirac delta function 
multiplying itself S(x) x 5(x). This product is not mathematically meaningful because its 
Fourier transform diverges: 

00 dk , 

— lxl — >oo, (8) 
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where we have used the convolution rule in Fourier transform and have noted that the 
Fourier transform of the Dirac delta function is 1. In contrast, the Hadamard method does 
allow multiplication for a wide class of distributions. For the example of S(x) x S(x), the 
Hadamard-neutralized Fourier transform of the product 



H °° dk 

H _„2? lxl = ' < 9 > 



yields 

5{x) x 6(x) = 0, (10) 

a mathematically meaningful (though somewhat counter-intuitive) result! 

In doing quantum field theory in configuration space, we multiply operator-valued dis- 
tributions of quantum fields. Or, in a slightly different interpretation, we multiply singular 
functions such as Feynman propagators. As we will see, the use of neutrix calculus allows 
one to put these products on a mathematically sound basis. Let us now generalize the 
above discussion for 1-dimensional Dirac delta functions to the case of (3 + l)-dimensional 
Feynman propagators 

, P A „ (11) 

(2tt)V + m 2 -ie K J 

where is the Hankel function of the second kind and we use the (+++-) metric. The 
Fourier transform of A + (x) x A + (x) (which appears in certain quantum loop calculations) 
is given by 

[d 4 xe- ip - x A + (x)A4x)=[-^-- \ — (13) 

J +K ' +K ' J (2tt) 4 k 2 + m 2 - ie (p - k) 2 + m 2 - ie V ; 



:D - 



4(2vr) 2 4(2tt)- 
where 

d 4 k 



f 1 dz]ag(l + £z(l-z)), (14) 
Jo \ m z J 



D 



A/ 

in J 



(k 2 + m 2 ) 2 
k 2 dk 2 



(15) 



'o (k 2 + m 2 ) 2, 

with the second expression of D obtained after a Wick rotation. But D is logarithmically 
divergent. Hence A + (x) x A + (x) is not mathematically well defined. Let us now see how the 



Hadamard-van der Corput method gives mathematical meaning to this product. Obviously, 
it is in the calculation of the logarithmically divergent D where we apply neutrix calculus. 
Introducing dimensionless variable q = k 2 /m 2 , we bring in to write D as 

where we have recalled that, for q — > oo, logq is negligible in the Hadamard neutrix H^. 
It follows that, in the neutralized version, A + (x) x A + (x) ~ S^(x)+ regular part (where 
S^(x) is the 4-dimensional Dirac delta function), a much more mathematically palatable 
object. 

As the first example in the application of neutrices to QED, let us consider the one- loop 
contribution to the electron's self energy 

= -ie 2 ( d " k ^[-7-Cp-fc)+m] 7 f ( 
[P) J (2vr) 4 [fc 2 + A 2 ][(p-A;) 2 + m 2 ]' 1 ' 

where m is the electron bare mass and we have given the photon a fictitious mass A to 
regularize infrared divergences. Expanding about 7 - p — —in, 

E(p)=A + B(yp + m)+K, (18) 

one finds (cf. results found in Ref.jsj]) 

A=-^-m(lD + -Y (19) 



2tt V2 4 

B=~^-4/:^ + ^, (20) 



47T V J±- x 2 

where a = e 2 /47r is the fine structure constant and D is given by Eq. and Eq. ()16|) 
in the pre-neutralized and neutralized forms respectively. We note that TZ, the last piece 
of S(p) in Eq. (fTSj) . is finite. Mass renormalization and wavefunction renormalization are 
given by m ren = m — A and ?/Ven = Z 2 i/j respectively with Z 2 =1 — 13. Now, since 
D = — 1 is finite, it is abundantly clear that the renormalizations are finite in the framework 
of neutrix calculus. There is no need for a separate discussion of the electron vertex function 
renormalization constant Z\ due to the Ward identity Z\ = Z 2 . 
The one-loop contribution to vacuum polarization is given by 

M*) = i(? / 7^k Tr U f AvT 7v , \, M ) ■ (21) 
J (27t) 4 \ 7-(p+f)+m 7-(p— f)+m/ 



A standard calculation P| shows that n M „ takes on the form 

n M „ = 5m\„ + (k\„ - k^K)Ii{k 2 ) , (22) 

where rj^ is the flat metric (+++-), 

5m 2 = — {m 2 D + D'), (23) 
2n 



and 

n(A; 2 ) = -— P + + — / dx^l-^log l + -x(l-i) , (24) 



with 



1 /■ <i 4 p 
m 2 J p 2 + m 

and .D given by Eq. Just as .D is rendered finite upon invoking neutrix calculus (see 

Eq. (PU)), so is -D': 

.D' = m 2 = 0, (26) 

Jo q + 1 V 7 

since both q and logg, for g — > oo, are negligible in if^. Thus neutrix calculus yields a finite 

— 1/2 

renormalization for both the photon mass and the photon wavefunction A^ en = Z 3 N 1 

1/2 l 

(and consequently also for charge e ren = Z s e) where Z% = 1 — 11(0). In electron-electron 
scattering by the exchange of a photon with energy-momentum k, vacuum polarization 
effects effectively replace e 2 by e 2 /(l — II(fc 2 )), i.e., 

e 2 

2 2 _ c 

6 6e// ~ i - n(fc 2 ) 

e 2 

ren 

" Z 3 (l-U(k 2 )) 

e 2 

^ren 

~ i - (n(fc 2 ) -n(o))' 

Eq. can be used, for k 2 3> m 2 , to show that 



(27) 



« e //(^ 2 ) = - - : 7 fc2 — • (28) 

1 3tt lQ S Vexp(5/3) 

Thus we have obtained the correct running of the coupling^] with energy-momentum in 
the framework of neutrices. In fact, the only effect of neutrix calculus, when applied to 
QED (and other renormalizable theories), is to convert infinite renormalizations (obtained 
without using neutrix calculus) to mathematically well-defined finite renormalizations. As 



6 



far as we can tell, all (finite) physically observable results of QED are recovered. In passing 
we mention that the use of neutralized integrals does not affect the results of axial triangle 
anomalies. 

As shown by the appearance of photon mass in the above discussion of vacuum polar- 
ization, the application of neutrix calculus to the energy-momentum cutoff regularization, 
though straightforward and natural, is ill suited for more complicated theories like those 
involving Yang-Mills fields. For those theories, one should use other more convenient regu- 
larization schemes. It is amusing to note that already in 1961 van der Corput suggested that, 
instead of finding the appropriate neutrices, one can continue analytically in any variable 
(presumably including the dimension of integrations) contained in the problem of tackling 
apparent divergences to calculate the coefficients of the corresponding asymptotic series. In 
hind sight, one recognizes that this was the approach taken by 't Hooft and Veltman who 
spearheaded the use of dimensional regularizationsj^|. Let us now explore using neutrix 
calculus in conjunction with the dimensional regularization scheme. In that case, negligible 
functions will include 1/e where e = 4 — n is the deviation of spacetime dimensions from 4. In 
the one-loop calculations for QED, the internal energy-momentum integration is now over n 
dimensions. The forms of S(p) and IT^ remain the same as given by Eqs. (fTHj) and but 
now with Sm 2 = 0. Using the approximation for the gamma function, T(e) = e _1 —7 + (9(e), 
where 7 ~ .577 is the Euler-Mascheroni constant, and the approximation f e ~ 1 + elog/, 
for e <C 1, one finds 

A = —[3(7 - log4vr) + 1] + — m / dx{l + x) \ogD , 
47T 2n Jo 

^ ot r , , am 2 f 1 , x(l — x 2 ) a f 1 , . , , „ , . 

B = — l+ T -log47T + / dx - 2 l / V + 77-/ dx(l-x)logD ,(29) 

4-7T f Jo m z x z + X z (l — x) 2tt Jo 

where D Q = m 2 x 2 + A 2 (l — x), and 

n(& 2 ) = 77- [7 — log 47r] + — ( dxx(l — x) log[m 2 + x(l — x)k 2 ], 
on ' n Jo 

I (X OL 

— = 1 - —[7 - log47r] - — logm 2 . (30) 
Z3 6n on 

By design, the generalized neutrix calculus renders all the renormalizations finite. Again, 
all physically measurable results of QED appear to be recovered. In this article we have 
explicitly considered QED to one-loop only. But we expect that higher-loop calculations 
can be handled in the same way according to neutrix calculus. It will be interesting to 
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see explicitly whether neutrix calculus, applied to higher-loop calculations, can provide new 
insights in the issue of overlapping divergences. 

In the framework of quantum field theory for the four fundamental forces, the divergence 
problem is particularly severe for quantum gravity. Using dimensional regularization, 'tHooft 
and Veltmanj7j found that pure gravity is one-loop renormalizable, but in the presence of a 
scalar field, renormalization was lost. For the latter case, they found that the counterterm 
evaluated on the mass shell is given by ~ e^^/gR 2 with R being the Ricci scalar. Similar 
results for the cases of Maxwell fields and Dirac fields etc (supplementing the Einstien field) 
were obtained p]. It is natural to inquire whether the application of neutrix calculus could 
improve the situation. The result is that now essentially the divergent e^ 1 factor is replaced 
by —7+ constant. 

It has not escaped our notice that neutrix calculus may ameliorate the hierarchy problem 
in particle physics. The hierachy problem is due to the fact that the Higgs scalar self-energies 
diverge quadratically, leading to a stability problem in the standard model of particle physics. 
But neutrix calculus treats quadratic divergences no different from logarithmic divergences, 
since both divergences belong to (the negligible functions of) the neutrix. Neutrix calculus 
may also ameliorate the cosmological constant problem in quantum gravity. The cosmologi- 
cal constant problem can be traced to the quartic divergences in zero-point fluctuations from 
all quantum fields. But again, neutrix calculus treats quartic divergences no different from 
logarithmic divergences. Indeed, for a theory of gravitation with a cosmological constant 
term, the cosmological constant receives at most a finite renormalization from the quantum 
loops in the framework of neutrix calculus. 

We conclude with a comment on what neutrix calculus means to the general question 
of renormalizability of a theory. We recall that a theory is renormalizable if, in loop cal- 
culations, the counterterms vanish or if they are proportional to terms in the original La- 
grangian (the usual renormalization through rescaling). It is still renormalizable if, to all 
loops, the counterterms are of a new form, but only a finite number of such terms exist. By 
this standard, neutrix calculus does not change the renormalizability of a theory, since it 
merely changes potentially infinite renormalizations to finite renormalizations. On the other 
hand, non-renormalizable terms, i.e., terms with positive superficial degree of divergence, are 
tolerated in neutralized quantum field theory. In a sense it is a pity that we have lost renor- 
malizability as a physical restrictive criterion in the choice of sensible theories. However, we 
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believe that this is actually not as big a loss as it may first appear. Quite likely, all realistic 
theories now in our possession are actually effective field theories. 0, 13 They appear to 
be renormalizable field theories because, at energies now accessible, or more correctly, at 
sufficiently low energies, all the non-renormalizable interactions are highly suppressed. By 
tolerating non-renormalizable terms, neutrix calculus has freed us from the past dogmatic 
and rigid requirement of renormalizability. (Having said that, given a choice between renor- 
malizable field theories and effective field theories, we still prefer the former to the latter 
because of the former's compactness and predictive power. But the point is that both types 
of theories can be accommodated in the framework of neutrix calculus.) Furthermore, if 
the application of neutrix calculus to loop calculations results in a term of a new form (like 
the Pauli term in QED) that is finite, then we have a prediction which, in principle, can be 
checked against experiments to confirm or invalidate the theory in question. For the latter 
case, we will have to modify the theory by including a term of that form in the Lagrangian, 
making the parameter associated with the new term an adjustable parameter rather than 
one that is predicted by the theory. This loss of predictive power is again not as big a loss 
as one may dread. 

Lastly we should emphasize that, for renormalizable theories as well as non-renormalizable 
theories (like quantum gravity?), neutrix calculus is a useful tool to the extent that it is 
relevant for asymptotic series and lessens the divergence of the theories. Based on our study 
so far, we tentatively conclude that neutrix calculus has banished infinities from quantum 
field theory, rendering perturbative quantum field theory mathematically meaningful. 
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